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Abstract 

Background: Neonatal cholestasis is a connnnon presentation of childhood 
liver diseases and can be a feature of various conditions including disorders of 
bile acid biogenesis and transport, various inborn errors of nnetabolisnn and 
perinatal infections. Sonne inherited nnetabolic diseases can be easily screened 
using biochennical assays, however nnany can only be accurately diagnosed by 
DNA sequencing. Fluorescent capillary Sanger sequencing (FS) is the gold 
standard nnethod used by clinical laboratories for genetic diagnosis of nnany 
inherited conditions; however, it does have linnitations. Recently nnicroarray 
resequencing (MR) has been introduced into research and clinical practice as 
an alternative nnethod for genetic diagnosis of heterogeneous conditions. In this 
report we connpared the accuracy of nnutation detection for MR with FS in a 
group of patients with 'low-nomnal' gannnna glutannyl transpeptidase (gGT) 
cholestasis without known nnolecular diagnoses. 

Methods: 29 patient DNA sannples were tested for nnutations in the ATP8B1 
and ABCB1 1 genes using both FS and MR. Other known causes of "low gGT 
cholestasis" such as ARC syndronne and bile acid biosynthesis disorders were 
excluded. 

Results: Mutations were identified in 13/29 sannples. In 3/29 sannples FS and 
MR gave discordant results: MR had a false positive rate of 3.4% and a false 
negative rate of 7%. 

Conclusions: The nnajor advantage of MR over FS is that nnultiple genes can 
be screened in one experinnent, allowing rapid and cost-effective diagnoses. 
However, we have dennonstrated that MR technology is linnited in sensitivity. 
We therefore reconnnnend that MR be used as an initial evaluation, with FS 
deployed when genetic and clinical or histopathological findings are discordant. 
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iij:j»7iii^ J Changes from Version 1 

In response to the referee comments, the discussion has been 
edited to acknowledge that several novel sequencing technologies 
have become available that are likely to or have already superseded 
microarray resequencing for clinical sequencing. These alterations 
help to put the work into context of recent advances in technology. 

See referee reports 



Introduction 

Neonatal cholestasis is characterised by persistent hyperbilirubi- 
naemia and has an incidence of around 1 in 2500 live births^ It can 
occur as a result of impaired bile acid biosynthesis, defective bile 
secretion by hepatocytes (due to membrane transporter defects), 
or extrahepatic obstruction to biliary flow. Patients present with 
jaundice, dark urine, pale acholic stools and hepatomegaly. Older 
infants may have pruritus, develop fat- soluble vitamin deficiencies 
and fail to thrived The aetiology of cholestasis is varied and 
includes multiple inherited causes, such as progressive familial in- 
trahepatic cholestasis (PFIC), Niemann-Pick disease type C and ar- 
throgryposis, renal dysfunction and cholestasis (ARC) syndrome^ 
Furthermore cholestasis can be subdivided into 'low-normal' and 
high gamma glutamyl transpeptidase (gGT) categories depending 
on the level of serum gGT activity. 

Accurate diagnosis is essential in order to inform decisions on treatment 
and management. As an example, familial intrahepatic cholestasis 1 
(FICl) protein deficiency caused by mutations in ATP8B1 (PFIC 
type 1) and bile salt export protein (BSEP) deficiency caused by 
mutations in ABCBll (PFIC type 2) have similar presentations, 
however liver transplantation is used more frequently in PFIC 
type 2 then type V\ Genetic testing allows accurate assessment of 
the genetic risk of cholestasis in families, in addition to providing 
definitive diagnoses. Clinical laboratories tend to use fluorescent 
capillary sequencing (FS; Sanger method), because it lends itself 
to automation and is highly sensitive. However, FS is limited by 
machine capacity and relatively high consumables costs. 

Microarray resequencing (MR) is an alternative sequencing method 
previously used to detect mutations in patients with intrahepatic 
cholestasis-. We have performed initial testing of a 300kb custom- 
designed microarray (Birmingham ReseqUencing Microarray, 
BRUMl) in patients with known mutations^. BRUMl includes 
probes to sequence two genes (ATP8B1 and ABCBll) involved in 
PFIC. This short research article describes further evaluation of the 
BRUMl microarray using samples simultaneously tested by FS and 
MR in patients without previously known mutations and assesses 
the advantages and disadvantages of MR in clinical laboratory use. 

Materials and methods 

DNA samples 

DNA samples from 29 children (15 females) who presented with 
intrahepatic cholestasis in the first year of life, low-normal range 
gGT activity without ARC syndrome or defects in bile acid bio- 
synthesis were used. ARC syndrome was excluded by clinical ex- 
amination^ and bile acid biosynthesis disorders were excluded by 
determination of urinary bile acid profile by electrospray ionisation 



mass spectrometry. Thus there was a high likelihood of detecting 
mutations in ATP8B1 and ABCBll in these samples. 

Polymerase chain reaction 

Polymerase chain reaction was used to amplify the coding exons 
of ATP8B1 and ABCBll plus the intron-exon boundaries prior to 
both FS and MR. The PCR reaction used 50ng of DNA, 0.5|iM each 
primer (Sigma), 0.5mM dNTPs (Bioline) and 2 units of Taq DNA 
polymerase (Biomix Red, Bioline UK) in Ix PCR buffer. 

Fluorescent Sanger sequencing 

5|al of PCR product was treated with ExoSap IT (GE Healthcare) 
for 30 mins at 37°C before being denatured at 85°C for 15 mins. 
Sequencing primers (at a final concentration of 0.5|iM) was added to 
Ijil of BigDye 3.1 and Ix BigDye buffer (Applied Biosystems) in a 
final volume of lOjil. The sequencing programme was thermal cy- 
cling for 34 cycles of denaturing at 95°C for 20 seconds, annealing 
50°C for 20 seconds and extension at 60°C for 4 minutes. Resulting 
products were precipitated using EDTA/Na acetate and ethanol and 
centrifuged at 4000 rpm for 30 minutes. Products were washed in 
70% ethanol before being air dryed, resuspended in 10|il of HiDi 
formamide (Applied Biosystems) and denatured. Resuspended 
products were analysed using a DNA analyzer 3730x1 (Applied 
Biosystems) and the resulting sequence traces were analysed using 
Sequence Analysis 5.2.2 (Applied Biosystems). 

Microarray resequencing 

PCR product quantitation, pooling, fragmentation, labelling and 
hybridisation and was performed according to the GeneChip 
Custom Resequencing Array Protocol V2.1 (Affymetrix). Arrays 
were washed and stained using a FS450 fluidics station before be- 
ing scanned with a GCS3000 7G scanner (Affymetrix). Intensity 
files were generated using AGCC (Command Console Vl.O, Affy- 
metrix) and processed in GSeq 4.1 (Affymetrix). Base calling as- 
sumed the diploid model and a quality score threshold of 2 (default 
settings were used for all other parameters). 

Results 

ATP8B1 and ABCBll were analysed using FS and MR. The se- 
quence variants detected in 13 of the samples (44.5% detection 
rate) are listed in Table 1 . Eight variants had previously been de- 
scribed in association with PFIC and ten variants were novel. Four 
of the variants were nonsense mutations (22%), two were splice site 
changes (11%) and the rest were missense (67%). No insertions or 
deletions were detected using the specifically designed probes. No 
mutations were detected by either technique in the remainder of the 
patients and therefore the cause of 'low-normal gGT' cholestasis 
remains undetermined in these cases. The negative results in these 
cases may be explained by whole exon deletions or duplications, 
intronic mutations affecting splicing, or promoter region mutations 
in ATP8B1 and ABCBll, as such mutations would not be detected 
by either strategy. In addition, it is possible that further genes are 
involved in the phenotype of neonatal cholestasis with low gGT. 

In sample 17, MR detected a variant not confirmed by FS, constitut- 
ing a false positive result. Two samples (6 and 22) had compound 
heterozygous changes in ATP8B1 detected by FS but not by MR. 
Both samples had two variants in close proximity; therefore we 
speculate that these are in cis and that each sequence variant has 
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Table 1.The sequence variants identified in this study. 



Sample Gene 


DNA change'' 


Protein change"" 


Zygosity 


Predicted effect 


Previously Microarray 
reported? resequencing 


Fluorescent 
Sanger sequencing 


2 




c.850G>C 


P.V284L 


Heterozygous 


Missense 


TUb 


Y 


Y 


4 




c.2178+1G>T 




Heterozygous 


Altered splicing 




Y 


Y 


6 


ATP8B 1 


c.1010T>G 


P.M337R 


Heterozygous 


Missense 


1 NUVUI 


N 


Y 


ATP8B 1 


c.1018A>G 


P.M340V 


Heterozygous 


Missense 


1 NUVWI 


N 


Y 


7 


ATP8B 1 


c.208G>A 


P.D70N 


Heterozygous 


Missense 


TUb 


Y 


Y 


9 


ABOB11 


c.2170G>A 


P.D724N 


Heterozygous 


Missense 


INUVUI 


Y 


Y 


12 


ABOB11 


c.2611-2A>T 




Heterozygous 


Altered splicing 


TUb 


Y 


Y 


17 


ATP8B 1 


c.1660G>A 


P.D554N 


Homozygous 


Missense 


T Wb 


Y 


Y 


ATP8B 1 


c.1564G>A^ 


p.D522N^ 


Heterozygous^ 


Missense 


1 NUVUI 


Y 


N^ 


18 


ABCB11 


c.290T>G 


P.L97* 


Homozygous 


Truncated protein 


Novel 


Y 


Y 


19 


ATP8B1 


c.3040C>T 


P.R1014* 


Heterozygous 


Missense 


Yes^ 


Y 


Y 




ATP8B1 


c.1014C>G 


P.N338K 


Heterozygous 


Missense 


Novel^ 


N 


Y 


22 


ATP8B1 


c.1018A>G 


P.M340V 


Heterozygous 


Missense 


Novel^ 


N 


Y 




ABCB11 


c.1636C>A 


P.Q546K 


Heterozygous 


Missense 


Novel^ 


Y 


Y 


25 


ABCB11 


c.499G>A 


P.A167T 


Heterozygous 


Missense 


Yes^2 


Y 


Y 


ABCB11 


c.3458G>A 


P.R1153H 


Heterozygous 


Missense 


Yes^° 


Y 


Y 


26 


ABCB11 


c.3484G>T 


P.E1162* 


Homozygous 


Truncated protein 


Novel 


Y 


Y 


29 


ABCB11 


c.483C>A 


p.cier 


Homozygous 


Truncated protein 


Novel 


Y 


Y 



^ - Variant was not confirmed by FS and is therefore a false positive finding. 
^ - DNA changes were experimentally determined by sequencing. 
° - Protein changes are predicted rather than experimentally determined. 
^ - Novel missense changes are of unclear pathogenicity. 

impaired the hybridisation of surrounding probes. Our experience 
of testing patients for mutations in these genes suggests that false 
negative results arising in this manner are likely to be uncommon. 
However, in this cohort they occurred in 7% of samples. 

Discussion 

MR is an attractive alternative to traditional FS for genetic testing 
in neonatal cholestasis. The BRUMl microarray used in this study 
has a much larger capacity than was utilised in this experiment, as 
it contains probes for 92 genes associated with inherited disorders. 
Analysis of the ability of this microarray to detect known gene 
mutations has shown a 97% mutation detection rate for base sub- 
stitutions^. The major advantage of MR is increased capacity, al- 
lowing sequencing of multiple genes in one experiment as quickly 
as one gene using FS. The main bottleneck is the requirement to 
quantify and to pool individual PCR products before hybridisation. 
This step is time-consuming and prone to error (due to pipetting 
volume variations and potential sample mix-ups). Therefore, au- 
tomation to minimise such errors would be useful, if not essential, 
for adoption of this method into clinical laboratories. Alternatively, 
use of long range PCR as the preparatory step for the MR protocol 
could reduce the number of reactions to be pooled, although it is 
more sensitive to poor DNA quality than standard PCR. Finally, 
microfluidics-based PCR systems, such as the 48.48 Access Array 
(Fluidigm Corporation, San Francisco, CA), might be combined 



with MR to avoid the quantification and pooling step altogether. 
This system allows the simultaneous but separate amplification of 
up to 48 PCR products for up to 48 samples, in nanolitre-sized re- 
actions, using a semi-automated process. The end product of the 
process is a pool of PCR products for each sample, and it is com- 
monly used for target enrichment for next generation sequencing 
experiments. 

Of the mutations recorded in the Human Gene Mutation Database^ 
in ATP8B1 and ABCBll, 24% and 16% were small insertions or 
deletions (indels) respectively, though their combined frequency 
in PFIC cases is unknown as most are private mutations. The major 
disadvantage of MR is it's insensitivity for the detection of indels^. 
Whilst known indels can be detected with proper microarray de- 
sign, novel indels will be missed. Larger insertions and deletions 
involving whole exons, of the type routinely detected by multiplex 
ligation-dependent probe amplification (MLPA; MRC-Holland, 
Netherlands), are not detected by either MR or FS. Another disadvan- 
tage, underscored by the results of this study, is that mis-called base 
substitutions (both false negative and false positive calls), are fre- 
quent and were found in approximately 10% of samples in this study. 

In summary, MR allows rapid and cost-effective genetic screen- 
ing in neonatal cholestasis and yields results relevant for pa- 
tient management. Many clinical laboratories are experienced in 
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oligoarray comparative genome hybridization (CGH) and thus 
have access to equipment required for MR, suggesting that MR 
could be implemented for a relatively small monetary invest- 
ment. In principle, MR could be applied to many clinical sce- 
narios involving heterogeneous conditions, especially if the 
mutations tend to be base substitutions. We recommend that MR 
methods be used for initial evaluation, with FS deployed when genetic 
and clinical or histopathologic findings are discordant. 

In recent years, various next generation sequencing (NGS) plat- 
forms have become available which allow massively parallel re- 
sequencing experiments to be performed. In particular 'benchtop' 
sequencers like the GS Junior (Roche), Personal Genome Machine 
(Life technologies) and MiSeq (Illumina) are attractive to clinical 
laboratories because they allow fast, cheap resequencing of multiple 
genes, and when the capacity is used optimally the costs are lower 
than MR. This study has identified a significant rate of discordant 
results obtained by MR when compared to Sanger sequencing, and 
consequently we predict that NGS will be the method of choice for 
clinical laboratory resequencing tests. 
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The authors compare, in a very thorough and detailed manner, the performance of IVIicroarray 
resequencing (IVIS) vs. Sanger sequencing (FS) for mutation detection in the ATP8B1 and ABCB11 
genes in cholestatic patients. They recommend the cautious use of microarray resequencing as an initial 
evaluation tool and propose to limit the use of Sanger sequencing to more problematic and discordant 
cases. 



However, given the non-negligible discordance observed by the authors between the two approaches 
(MR had a false positive rate of 3.4% and a false negative rate of 7%), the definition of the cases that 
should be analyzed by FS is not entirely clear. Moreover, a number of other resequencing procedures that 
can be applied to the simultaneous analysis of multiple genes are now available. Some of them are very 
competitive in terms of efficacy and cost with MS and it is quite likely that in the near future they will be 
more widely used with respect to MS. Ideally, the authors should have discussed about this issue. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard, however I have significant reservations, as outlined 
above. 



Competing Interests: No competing interests were disclosed. 



Struan Grant 

Center for Applied Genomics, Children's Hospital of Philadelphia Research Institute, University of 
Pennsylvania, Philadelphia, PA, USA 



Approved: 14 February 2013 
Referee Report: 14 February 2013 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 

Competing Interests: No competing interests were disclosed. 



Page 7 of 7 



